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Abstract 



Bjorken scaling is violetd. At large x-values (0.7 ^ x < 1) the violation is 
mainly attributed to oc l/Q"^ (higher-twist) corrections. We discuss how to 



> 

\ incorporate such corrections by using a new scaling variable x{x,Q'^) which 

O 



accounts for the non-perturbative effects due to the confining parton interac- 



\ tion in the final state. The approach can accomodate also the remaining QCD 

Qh! logarithmic corrections and the behaviour of the F2{x^ Q"^) structure functions 



is reproduced, in a quite natural way, within a wide kinematical range. 

The existing data for hadron structure functions, F2(x,(5^), show considerable Q^- 
dependence, which is mainly attributed to the QCD logarithmic corrections to Bjorken 
scaling. However, at x ^ 1 the scaling violations are dominated by power corrections 
oc (higher twist and target mass effects): 

F2(x,Q2)=Fr(x,Q2) + ^ + ---, (1) 

where F2^{x,Q'^) = F2{x,Q^ ^ |i?(a;)|) and the remaining Q^-dependence in F2^{x,Q'^) is 
to be attributed to QCD logarithmic corrections only. 

Power corrections can be incorporated in the first term of Eq. ([I|) by using a different 
scaling variable, 



X = Q^) = X + 



b{x) 



+ ■■■, 



(2) 



so that 



tas 



(3) 



The coefficient B, which determines the value of the power correction in Eq. is thus 
related to the structure function by B{x) = b{x) dF2'^{x, Q'^)/dx. 

Actually, an analysis of data in terms of an appropriate scaling variable appears to be 
more convenient, than the direct evaluation of power corrections. For instance, it is common 
to use for an analysis of data the Nachtmann variable x = ^ |Q, 



which is expected to account for the effects of the target mass (M). 

Besides the target mass effects, there are important nonperturbative effects from the 
confining interaction of the partons in the final state. Indeed, the partons are never free, so 
that the system possesses a discrete spectrum in the final state. Although in the Bjorken 
limit the struck quark can be considered a free particle, the discreteness of the spectrum 
manifests itself in power corrections to asymptotic structure functions 0. One can anticipate 
that these corrections are significant in particular at large x, where lower-lying excitations 
should play an important role. 

A general analysis performed in the framework of Bethe-Salpeter equation shows that in 
the case of a local confining final state interaction the higher twist and target mass effects 
can be effectively accounted for by taken the struck quark with the same off-shell mass 
before and after the virtual photon absorption 0. As a result, the Bjorken scaling variable 
X is replaced by a new scaling variable x = x{x, Q"^), which is the light-cone fraction of the 
off-shell struck quark. Explicitly, 



(4) 



X = 



(5) 
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where M is the target mass and iris is the invariant mass of spectator partons (quarks and 
gluons). For ^ oo or for a; — the variable x coincides with the Nachtmann variable ^, 
Eq. (^. However, at finite these variables are quite different. 

It follows from Eq. (^) that x depends on the invariant spectator mass, m^. The latter 
can be considered a function of the external parameters only. In the limit x ^ 1 (elastic 
scattering) no gluons are emitted, and thus tuq, which is the mass of a two-quark 

system (diquark). When x < 1, the spectator mass increases due to gluon emission. For 
X close to 1, nis can be approximated 



ml + Cil-x), (6) 



where the coefficient C ~ (GeV)^. In the following we regard it as a phenomenological 
parameter, determined from the data. 

Let us consider the nucleon structure functions in the region of large x, where the power 
corrections to the scaling are dominant. At present, the only available large- a; data for 
proton and deuteron structure functions, F2 {x , Q"^) , F2 {x , Q'^) , are the SLAG data 0-0, 
taken at moderate values of momentum transfer, ^ 30 (GeV/c)^. (The nucleon structure 
functions for higher values of momentum transfer {Q^ ^ 250 (GeV/c)^) are extracted from 
BGMDS H and NMG g] data, yet only for x < 0.75). The SLAG data for the proton and 
deuteron structure functions for x > 0.7 and 5 ^ ^ 30 (GeV/c)^ are shown in Fig. 1 as a 
function of x. Also shown is the value of F^{x, Q^) and F2{x, Q"^) for Q'^= 230 (GeV/c)^ and 
X = 0.75 taken from the BGDMS data P). The data points close to the region of resonances 
were excluded by a requirement that the invariant mass of the final state (M + z/)^ — 
is greater than (M + A)^, where A = 300 MeV. In addition, we excluded the data points 
with X > 0.9 from the deuteron structure only (Fig. lb). The reason is that the deuteron 
structure function can not be represented as an average of the proton and neutron structure 
functions for x ^ 0.9. Indeed, the calculations of Melnitchouk et al. |T0| show that the ratio 



2F2/{F^ + F2 ) is about 1.13 for x = 0.9 and Q^=5 (GeV/c)^, and it rapidly increases for 
X > 0.9. However, for x < 0.85, this ratio is within 5% of unity |jll . 



One finds from Fig.l that tlie structure functions show no scaling in the Bjorken variable 
X. Also, very poor scaling is obtained when the data are plotted as a function of the 
Nachtmann variable ^, Eq. However, the situation is different if we display the same 

data as a function of the variable x, Eq. (^). It appears that the scaling in the a;- variable is 
strongly dependent on the value of diquark mass, mo, Eq. (^, but is much less sensitive to 
variation of the coefficient C. For instance. Fig. 2 shows the data as a function of x for mo = 
600 MeV and C = 3 (GeV)^), i.e. by considering the spectator as build up from constituent 
quarks. The data display very poor scaling, although it is slightly better than that shown 
in Fig. 1. The scaling deteriorates even more when mo > 600 MeV. On the other hand, 
the scaling is very good both for the proton and deuteron data, by taking mo = 0, i.e. by 
considering the spectator build up by current quarks 0. The results are shown in Fig. 3, 
where the data are plotted as a function of x for mo = and C = 3 (GeV)^. Note, that the 
high-Q^ data points from BCDMS data P] are very close to the SLAG data points, taken 
at much lower values of Q^. Also note that x ^ 1 when a; 1 for mo = 0, Eq. but 
x{x, Q"^) < X for X < 1. As a result the data, plotted as a function of x-variable are extended 
in a wider region, than the same data plotted as a function of the a;-variable (cf. Figs. 1 - 
3). 

Now by using the scaling variable x, Eq. (H) for mo = we are going to analyze the proton 
structure function for smaller values of x, where both power and logarithmic corrections to 
the Bjorken scaling play an important role. In this region (0.35 < x < 0.75) the existing 
BGDMS [|] and NMG data are extended up to much larger values of momentum transfer 
than the previously considered high-x SLAG data. It allows us to check our predictions in 
a wide range. 

QGD (logarithmic) evolution effects on F2 are taken into account at Next-to-Leading 
Order (NLO) ||12[ evolving back, in Q^, the structure functions starting from an asymptotic 
value of momentum transfer where the condition F2{x,Q'^) ~ F2^{x,Q'^) (cf. Eq. ([1|)) is 
fulfilled (in the present case we choose = 230 (GeV/c)^, which is the highest value 
of the momentum transfer in the BGDMS data |^). At that value of the functional 
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form of F2 is parametrized following the prescriptions of the recent NMC fit and valid 
for a wide range of x (0.006 < x < 0.9). For the x-region we are interested in, namely 
X > 0.35, one can assume that the valence contribution to F2 dominates and one can consider 
evolution of the nonsinglet (NS) components only. The accuracy of such an assumption can 



be deduced from the recent CTEQ parametrization [III] of the parton distributions. The 
ratio r = F^^{x,Q^)/F2{x,Q^) is r ^ 0.96 for x = 0.35 and r > 0.995 for x > 0.55 and 
> 5 (GeV/c)^. In fact, since the only assumptions we are making is due to the fact that 
the (small) singlet component evolves in in a different way with respect the dominant 
nonsinglet part, the inaccuracy is less than 4% for x = 0.35 and less than 0.5% for larger 



X. For the actual calculations we use the NLO procedure developed in ref. |T^. Under the 
renormalization group equation (RGE) the moments of the nonsinglet components of the 
nucleon structure function, {F^^{Q'^))n = Jq F^'^(x, Q^), evolve according to 

{F2^'m)n = {F2''HQl))n ^ JZ,., Z (7) 



where 



l + ^RNs\asm 



l,n n 

tdNS 'Ins Pi o,n . ^i,n /q\ 



and we used the expansion 

as{Q^) 1 /3i lnln(QVA2 



(9) 



47r /3oln(gVA2) Pl[\n{Qy A^)]^ ' 

/3o = 11 — 2//3, (3i = 102 — 38 //3 and / is the number of active flavours. The values of 
all the other parameters of the calculation (like the anomalous dimensions 7)^5, 7]^^ and 
the Wilson coefficient 6*2^'^) have been compiled in refs. [|T2],|IB|. The form (|^ guarantees 
complete symmetry for the evolution from to > Qq and back. In our case Qq = 
230 (GeV/c)2 and 5 < < 230 (GeV/c)^. The NLO factorization scheme imphcitly 
selected in Eq. (^ is the so called DIS scheme where the structure function assumes the form 
Fiix, Q') = X {q{x, Q') + q{x, Q')). 
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The procedure of perturbative QCD evolution previously sketched implies that the power 
corrections are already extracted from the structure functions, so that the evolution has to 
be actually applied to Eqs. (|l|)-(@) 

F^ix, Q') = F,{x,Q')-^----c^F, {x-\x, Q'), Q') , (10) 

where we used the scaling variable Q^) = x{x,Q'^), Eq. (^) for mo = and C = 3 
(GeV)^ that effectively incorporates the power corrections. It implies that the evolution of 
structure functions feacA; from = 230 (GeV/c)^, should be calculated with a shifted value 
of Bjorken variable, namely 

Ai(x, = F2 {x-\x, Q^), Q^) - F2 {x-\x, Q^), 230) . (11) 

The influence of this shift would be quite important for low and moderated values of Q^. 

The corresponding variation of structure functions with for fixed x due to power 
corrections can be evaluated as 

A^ix, Q2) = F2 {x{x, Q^), 230) - F2 {x{x, 230), 230) . (12) 

Finally the Q^-dependence of structure functions due to logarithmic and power corrections 
to Bjorken scaling is given by 

F2{x, Q2) = 230) + Ai(a;, Q^) + A2(x, Q^). (13) 

The results are shown in Figs. 4 and 5 for proton and deuteron structure functions respec- 
tively. The data points are from SLAG and BGDMS data bins . The dotted lines show 
the Q^- dependence of the structure functions due to power corrections only, Eq. (p!2D. The 
total Q^-dependence of structure functions due to the power and the logarithmic NLO cor- 
rections, Eq. (p!3D, is shown by the dashed and continuous lines for A = 100 MeV and A = 200 
MeV respectively. One finds from Figs. 4 and 5 that Eq. (|T3|) reproduces the experimental 
data in a large Q^-range for both values of A, although the agreement is slightly better for 
A = 100 MeV. In addition, since the results are strongly dependent on the spectator mass 
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(^), it is remarkable that the same parameters mo = and C = 3 (GeV)^ do reproduce the 
Q^-behavior of the structure functions both for large and moderate x-values. 

In conclusion we have presented a detailed investigation of the F2 experimental data, both 
for proton and deuteron, analysed by means of the scaling variable x recently proposed in 
ref. 1^. Such variable includes the non-perturbative effects due to the confining interactions 
of the partons in the final state and contains, in an effective way, higher-twist corrections. 
These contributions show up at large x and at low and moderate momentum transfer Q^. 
The scaling behavior of the experimental data is improved in a wide range of x and Q^. The 
additional inclusion of the QCD radiative logarithmic corrections allows us the investigation 
of the structure function in an even larger kinematical range with the extra bonus of a rather 
precise identification of the spectator diquark mass. 
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FIGURE CAPTIONS 



Fig 1. The SLAC data {5 < < 30 (GeY/cf) for proton (a) and deuteron (b), 

are shown as a function of the Bjorken variablle x. Three high-statistics data sets |^ for 
Q2 ~5.7^ 7.6, and 9.5 (GeV/c)^ are marked by "+", "x", and "#" respectively. The point 
at g2 = 230 (GeV/c)2 and x = 0.75 is from ref. ||]. 

Fig 2. The data of Fig. 1 are shown as function of the x{x, Q^) — the scahng variable of 
Eq. (H) — assuming mo = 600 MeV and C = 3(GeV)^ for the spectator mass m^, Eq. (|]). 

Fig 3. As in Fig. 2 assuming mg = and C = 3 (GeV)^. 

Fig 4. The proton structure function -^^(x, Q^) is shown as a function of at different 
x-values. The dotted lines include power corrections only. They are evaluated according to 
Eq. ( [T^ ) and the scaling variable x of Eqs. (^), (|^) with rriQ = and C = 3 (GeV)^. The 
additional QCD logarithmic corrections evaluated at NLO according to the procedure of 
Eq. (]TT|), ( prSD for different A scales are shown by the dashed (A = 100 MeV) and continuous 
lines (A = 200 MeV). 

Fig 5. As in Fig. 4 for the deuteron structure function F2{x,Q'^). 
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